Submitted to Physical Review Letters 



Many— Body Corrections to Charged— Current Neutrino 
Absorption Rates in Nuclear Matter 



Adam Burrows 

Department of Astronomy, The University of Arizona, Tucson, AZ 85721 
e-mail: burrows@as.arizona.edu 



R. F. Sawyer 

Department of Physics, The University of California, Santa Barbara, CA 93106 
e-mail: sawyer@sarek.physics.ucsb. edu 



Abstract 



Including nucleon-nucleon correlations due to both Fermi statistics and 
nuclear forces, we have developed a general formalism for calculating the 
charged-current neutrino-nucleon absorption rates in nuclear matter. We 
find that at one half nuclear density many-body effects alone suppress the 
rates by a factor of two and that the suppression factors increase to ~5 at 
4 x 10 14 g cm -3 . The associated increase in the neutrino-matter mean- 
free-paths parallels that found for neutral-current interactions and opens up 
interesting possibilities in the context of the delayed supernova mechanism 
and protoneutron star cooling. 
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I. INTRODUCTION 



The neutrino absorption and scattering opacities in the post-shock core of a supernova, in 
which nuclei are largely disintegrated into nucleons, determine the duration, spectrum, and 
flavor distribution of the emerging neutrino pulse. It has been known for some time that the 
interactions among nucleons in the denser regions can change these opacities significantly, 
but to date there has been no comprehensive treatment given in the literature and present 
calculations of the complete supernova process do not include the effects of interactions on 
the opacities. 

The neutrino-matter interaction rates can be related to the space- and time-dependent 
correlations among the set of density operators for the separate nuclear constituents (to find 
the Gamow-Teller parts we must consider separate spin-up and spin-down densities). In 
the case of neutral-current interactions [0J, there is an instructive limit, which also provides 
an estimate of the effects, in which the combined limits of large nucleon mass and small 
neutrino energy allow the use of long-wavelength limits of equal-time correlation functions, 
in turn expressible in terms of the second derivatives of an energy density functional with 
respect to various densities. This approach is the direct multichannel generalization of the 
familiar results for light scattering from the thermal density fluctuations in a fluid, where 
it is the compressibility that determines the long- wavelength opacity, and it was used in 
references [0 and || to find significant reductions of neutral-current opacity in certain 
regions. In Burrows & Sawyer [[I]] , an approach based on ring graphs was used to encompass 
these results and to extend them to domains in which the equal-time and long-wavelength 
limits are not clearly applicable. 

The use of the equal-time and long-wavelength limits to express correlation functions in 
terms of static susceptibilities cannot be extended to the charged-current interactions when 
there is a large chemical potential difference between protons and neutrons. Furthermore, 
there do not exist in the present literature systematic estimates of the effects of interactions 
on the charged-current opacities for electron neutrinos. In the present work, we give a 
theoretical framework for addressing these opacities, based on summing ring graphs, together 
with the results of calculations with input parameters taken from the current phenomenology 
of nuclear matter. 

II. MANY BODY FORMALISM FOR CHARGED CURRENT RATES 

The charged-current interactions of the nuclear medium are determined by the statistical 
averages involving the nucleonic charged-current operator, = T/yy^l — gAl5)4>n, and its 
Hermitean conjugate, 

W, u {q^) = -iZ- 1 j dW^V^2V[ e -^-^ (1) 

where Z is the partition function and ft is the chemical potential difference, fi = fi n — fi p . 
The vector, (q, u), is the momentum-energy transfer from the leptons to the medium. The 
transition rate, differential in angle and u, for the reaction v e + n — > e~ + p is given in terms 
of W au as 
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^T.l) = (^)-^ J E \- f l ) - f ^-^ A^ M ^W M (2) 
dudcos{0) 1 — exp|— p(u + fi)\ 

where 

A"" = (4^^)-^^ (1 - 75)7^2 7^(1 - 75)] , (3) 

q = [Ef + (Ei — to) 2 — 2Ei(Ei — cu) cosO] 1 ^ 2 , and E\ in the incident neutrino energy. To 
verify the factors in (0) replace the commutator in (|1|) by the unordered product, the 8(t) 
function by 1/2 in order to capture the imaginary part, and the statistical factor in the 
denominator by unity. Then, the introduction of a complete set of states between the 
current operators gives the inclusive v e —>■ e~ rate. Recalling that the Heisenberg picture 
for the density operators is defined with respect to the operator H — E/XjATj, we note that 
the e 1 ^ factor in ([]]) will be canceled before the time integral is performed. The other term 
in the commutator gives the same answer, but with an additional factor, — exp[— j3(u + jl)}. 
The statistical factor in the denominator of (0) cancels when the two terms are combined. 

Taking the nucleons as nonrelativistic, which is allowable up to about twice nuclear 
density, we note that the components Wio and Woi vanish, and that the space tensor, Wij, 
derives from the axial-vector current and can be written in the form, 

Wij(u, cosO) = g\WA(uJ, cos 9)5ij + g\Wr(uj, cos 9)qiqj . (4) 

We define Woo — Wy, as it comes entirely from the vector current. Combining with (^) and 
calculating the lepton trace, we now have 

W^A^ = 2(1 + cos6)W v + 2(3 - cos6)g 2 A W A + 2uj 2 (1 + cos6)g 2 A W T . (5) 

We calculate the three functions Wy, W A , and Wt in the ring approximation. For these 
graphs, a unit of charge is passed in each nuclear interaction in the expansion of the charged- 
current correlation function of (0). If the potential between nucleons a and b conserves 
isotopic spin, this means that only the terms containing the isotopic operator 
enter. We assume a potential in momentum space, 

V a ' b (q,uj) = r a ■ r b [ Vl (q) + a a ■ a b v 2 (q) + a a ■ qa b ■ qv 3 (q,u)\ . (6) 

We have inserted in the V3 term a dependence on energy transfer, to, in order to acco- 
modate the one-pion exchange force, 

v 3 (q,co) = -f 2 m- 2 (q 2 + ml-u J 2 )- 1 , (7) 

with f«l. 

For the other two potentials, we follow the development of Jl] in taking zero-range forms 
fitted to the Landau parameters of reference j|] and obtain t> 2 = 3.4m/m* x 10~ 5 (MeV)~ 2 
and v\ = 1.88m/m* x 10 _5 (MeV)~ 2 , where m* is the effective nucleon mass in the medium. 
We are assuming isospin invariant forms, despite the fact that our application will be to 
unsymmetric matter. This appears to us to be the state of the art. It is not totally satis- 
factory, in view of the fact that the v 1 and v 2 terms are phenomenological forms that obtain 
in nuclear matter rather than forms based on elementary meson exchanges. Note that in 
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the case of non-symmetric matter isospin symmetry is still broken through the polariza- 
tion functions. We define the nucleon-charge-raising polarization part, U(q,u), in parallel 
with (Jl]), in terms of the retarded commutator of the density n c = ippjjn and its Hermitean 
conjugate, 

U{q,u) = -1Z- 1 J d A x e~ i(l - x e l{uJ+fi)t Tr[e-^ H - E ^ N ^[n c (x, t), 4(0, 0)]]6(t) . (8) 

Note that the function that gives the vector part of the rate, Wy(q,u), is exactly given by 
II. The ring approximation |j| for II is now 

WvM = U M = T -^^ ry (9) 

where II^ ) is the polarization function in the absence of interactions. The factor of 2 
multiplying the potential comes from the isospin operator in the potential ([]). The spin 
dependent parts of the potential do not contribute. 

For the axial contribution, the operative non-relativistic forms of the current are the 
operators n l c = ip^ipn, and the analogue to (§) will be a tensor with indices In the 
absence of interactions, this tensor is given by 5 4j n< '. Thus, summing the axial chain to 
get WA{q, oj) we obtain exactly the same structure as (|8|) with v 2 replacing v\, but with the 
same function 11^, 

w = i- 2 Sn*Lr (10) 

The potential does not enter fllOD, even though it couples to the axial-vector terms, since 
any ring graph chain in which at least one participates becomes a contribution to Wt- 
Elementary combinatorics for the tensor chain gives, 

w T , q u) = [n<%.^M 

T[q ' ' i-2H!)+W!,»)]n |, i(*,«) 

The requisite polarization function is given by 

n( 0) ( } = _ 2 /-^P /(|pUn)-/(| P + q|,/x P ) (12) 
J (2ny oj + e p - e p+q + ir] 

where the functions / are the nuclear Fermi occupation functions for the indicated momenta 
and chemical potentials. We conceptually extend fll2"D to include the average potentials that 
the nucleons experience in the medium, v pM , by making the replacements (p) 2 /(2m) — > 
(p) 2 / (2m) + v n and (p + q) 2 /(2m) — > (p + q) 2 /(2m) + v p , both in the denominator and in 
the distribution functions in the numerator. However, we suppose that we are starting with 
a table of densities and temperatures (T) from an equation of state that already takes into 
account the potentials v PjTl . If we utilize Fermi distributions in which the chemical potentials 
are derived from the input densities using the free-particle relations, these average potential 
corrections are automatically included and the parameters v p>n do not appear explicitly in 
the formalism. To do the computations, we use the following form for II^ ^ : 
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(15) 



+ (a; -> /i n -> /i p ) . (16) 



The imaginary part of the polarization (|15|) can be obtained by direct integration, and is 
the same as that given in references |1]||. Given this, the full function ( |T3"D can be verified 
by checking the analytic properties in the uj plane. 

We look at the results of medium interactions for two sets of conditions typical of the 
dense interior of a protoneutron star, but at different times; first for an early time (t < 5 
seconds), when the lepton number is large and second at a later time (t > 10 seconds), after 
which the lepton excess has largely radiated away M. 



III. LEPTON RICH ERA 

The densities of the various species are such that at zero temperature the inequality 
2pp > — pp | holds. This inequality allows the single nucleon process to proceed at 
full strength. That is to say, at low temperatures the function ImlT '' (q, ui) flT5] ) is large in 
the region of (q, oj) defined by the two conditions: 1) the neutrino energy, E%, is near the 
neutrino Fermi surface; 2) the electron energy, E\ — u, is near the electron Fermi surface. 
These conditions come from the occupation factors that enter the expression for the total 
rate. 

In this region, we use (Q), fllQ|) , and (|TID to calculate the modification factors in the 



medium. We have estimated the contribution of the tensor term (|TTD relative to the two other 
terms and conclude that it is less than 10% of the total in the cores of protoneutron stars 
and supernovae. For the vector and axial-vector terms, we compute suppression factors, 
5*^4 and Sy, defined as the ratio of the rates calculated with the nuclear interactions to 
those calculated without. This is done by substituting @ and ( |10"D into (H), multiplying 



(g) by the neutrino occupation function, and integrating over neutrino energies. In Table 
1, we give results using a post-bounce supernova profile taken from Burrows, Hayes, & 
Fryxell [§]. As seen in the table, the Gamow-Teller suppression factors are larger than the 
Fermi suppression factors. Furthermore, the degree of many-body suppression increases 
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with density and decreases with temperature. Importantly, the magnitude of the effect 
above 10 14 g cm -3 is large, ranging from a factor of 2 at 10 14 g cm -3 to a factor of ~5 near 
4 x 10 14 g cm -3 . Correspondingly, the neutrino-matter absorption cross sections decrease 
with density. Since it has recently been shown |l| that the neutral-current scattering rates at 
high density are also reduced, we conclude that post-bounce supernova cores are significantly 
more transparent than previously believed. This enhanced transparency should translate at 
late times (> 500 — 2500 ms after bounce) into higher neutrino luminosities, that thereby 
may be more efficient at reenergizing or powering a stalled supernova shock [|[]. 

Integrating equation (^) over cos6, the distribution of the energy transfer, u, to and from 
the nucleons due to the process u e + n — > e~ + p in the lepton-rich era can be derived and 
is depicted in Figure 1 for a variety of densities, from 10 13 g cm -3 to 10 15 g cm -3 . For these 
curves, the temperature is 5 MeV, the electron fraction, Y e , is 0.26, the incident electron 
neutrinos are on their Fermi surfaces, and beta equilibrium is assumed. The highest curve 
on the right (that for p = 10 13 g cm -3 ) ignores many-body effects, though it incorporates 
the full kinematics, and is included for comparison. As expected, the peak of the energy 
transfer is generally near —ft (given in the figure caption), since the electron blocking factor 
in (0) puts the electrons on the electron Fermi surface and beta equilibrium requires that 
\i e = ft + fi Ue . (Note that ft increases with density.) There is a modest spread in to around 
the peak with approximately a gaussian distribution. The width of this distribution scales 
with the temperature. Figure 1 demonstrates what Table 1 also reflects that the total cross 
sections, the integrals under the unnormalized curves, are decreasing functions of density. 



IV. LEPTON POOR ERA 



As discussed in || and fllU| , as the trapped electron lepton number decays, we reach a 
configuration in which the neutrino absorption process discussed above dies almost com- 
pletely for low temperatures. At the end of deleptonization, we have pp = [f|j. The 
neutrinos then have thermal energies and we find that the proton fraction has decreased 
to the point that 2p^p << pffi. In this case we cannot conserve momentum for the three 
degenerate species, e~,p,n, when the momentum of the neutrino is small and when we 
stay near the Fermi surface for the three other species. At low temperatures, the function 
Imn(°)(g,cj) of (0) is now exponentially small in the region of (q, to) defined by the leptonic 
occupancies. Thus, the ring graphs, as defined in the previous section, give a negligible rate 
for the charged-current processes at low lepton number and temperature. 

The other mechanisms that have been used to estimate the rates depend on a spectator 
nucleon to transfer the necessary momentum, either through a potential or through an 
assumed correlation P-pJjj. Translated to graphs, these mechanisms involve the estimation 
of proper graphs for the polarization parts, where a proper graph is defined as a graph that 
cannot be cut into two disjoint parts by severing a single potential line. However, there are 
"ring corrections" to such graphs, in which the initial or final current vertices attach to a 
ring chain that then attaches to the proper polarization graphs. These ring corrections then 
have a big suppressive effect on the primary mechanism for momentum transfer, as one can 
see from the following argument. 

We consider a proper polarization graph in which a nucleon-nucleon interaction has 
intervened to allow momentum conservation near the Fermi surfaces. That is to say, by 
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([10|), we have a term for the vector contribution, Wy, to the rate formula (|5|) that has a 
substantial imaginary part in the kinematically-allowed region. We call this term, Wy . We 
now take the sum of this term and the lowest order term, U^°\ as the proper polarization 
part from which to construct the ring sum, obtaining, 



ImW 



V 



^y( q ,u) « - - — - V.^)' ,(,), - • (17) 



, - 2v 1 [UW(q,cu) +ReW v I> + ilmWy 

In the estimates that follow, we omit the real part of Wy from the denominator; it is easy 
to verify from (|T6| ) that, in contrast to the lowest-order imaginary part, the real part is not 
suppressed in the region of (q, uj) that is important in the reaction. 

Similar considerations hold for the axial- vector part. In the lepton-poor era, we can 
drop the second term on the RHS of (|5|), since uj will be of the order of T, rather than of 
order 100-200 MeV, as it can be for electron neutrinos in the trapped neutrino era fll% 
Then, we define a proper Wa and a contribution W A as i n the above and write 



\1 - 2v 8 [UW(q,u) + ilmW\ 1 \q,uj) 



For the factors ImW 7 ^ we take the minimal form consistent with avoiding a singularity 
in (Q) at uj — —ft: 

lmW$ A = c$ A (uj + fi). (19) 

This form meets the requirement for detailed balance that Imll be odd under the replace- 
ments, uj — > — uj, (1 — > —ft. We take the parameters c- 1 ' to be sufficiently small for the 
integrated suppression factor to be independent of the [c^] 2 that enter through the imagi- 
nary parts in the denominator function. In Table 2 we give, for the case of matter at nuclear 
density, the separate suppression factors for the vector and axial-vector rates generated by 
the uncorrected terms (|T9"D, under deleptonized conditions, //„ = 0. We note that the reduc- 
tions are substantial, and conclude that the ring corrections should be added to any model 
that is used for the neutrino opacities during this era. 



V. CONCLUSIONS 

We have developed a new formalism for incorporating the effects of many-body correla- 
tions on the charged-current rates of neutrino-matter interactions. This formalism reveals 
that these rates are considerably suppressed in the densest regions of protoneutron stars and 
supernova cores. Assuming that the nucleons are non-relativistic, our formalism incorpo- 
rates the full kinematics of the interaction, Pauli blocking by final-state nucleons (protons), 
and correlation due to nucleon-nucleon interactions. 

We have employed the ring approximation (RPA) and assumed the near-validity of 
Fermi Liquid Theory. It would desirable to include ladder diagrams and to perform the 
calculations in the context of a better numerical method for solution of the nuclear equation 
of state (EOS), since the solution of the EOS is intimately related to the derivation of the 
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scattering/absorption rates. However, those who perform detailed nuclear EOS calculations 
and address many-body correlations in nuclear matter do not as yet provide the requisite 
spin and density structure functions, even for the static case. 

These results for charged currents, when combined with the results from Burrows & 
Sawyer |I[ for neutral currents, strongly suggest that energy and lepton number will leak 
from supernova cores at a rate that is higher than heretofore estimated. This implies that 
the neutrino luminosities during the epoch after bounce for which the inner core is the 
major energy source (> 0.5 — 1.5s) will be enhanced, perhaps by as much as 50% |IJ. The 
consequences of this increased transparency for the neutrino-driven supernova explosion 
mechanism || may be interesting, but have yet to be clarified. 
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TABLES 



TABLE I. The total suppression factors (<S,4,v) f° r the process v e + n — > e _ + p, for a profile 
in an early post-bounce model generated by Burrows, Hayes, & Fryxell (1995). The suppression 
factors for the vector and the axial-vector terms are shown separately. These suppresion factors are 
derived by multiplying the rate by the neutrino occupation function and integrating over neutrino 
energies. A nucleon effective mass of 0.75 x m n is assumed. 



P (g cm" 3 ) Y„ T (MeV) Sa Sv_ 

3.94 x 10 14 0.077 5 0.289 0.140 0.269 

3.68 x 10 14 0.078 5 0.294 0.144 0.275 

3.08 x 10 14 0.077 5 0.297 0.157 0.291 

1.65 x 10 14 0.064 10 0.275 0.228 0.381 

2.66 x 10 13 0.01 15 0.282 0.670 0.775 
1.40 x 10 13 0.067 15 0.258 0.790 0.840 



TABLE II. The total axial and vector suppression factors (<S4,y) for the transformation v e — > e~ 
during the lepton-poor era, for a density of 2.5 x 10 14 g cm~ 3 and a neutrino chemical potential 
of zero. These suppression factors are derived by multiplying the rate by the neutrino occupation 
function and integrating over neutrino energies. A nucleon effective mass of 0.75 x m n is assumed. 
See the text for details. 



T (MeV) Y. Sa Sv_ 

3 0.012 0.34 0.53 

5 0.013 0.20 0.38 

7 0.016 0.19 0.38 

9 0.018 0.20 0.38 
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creating a proper polarization part that has a significant imaginary part in the region 
of interest, for either the vector or axial-vector current parts. 
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FIGURES 
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FIG. 1. The singly-differential cross section in arbitrary units versus the energy transfer, lu, to 
the nucleons due to the process v e + n — > e~ + p in the lepton-rich era. u> is in MeV. We have 
integrated Eq. (2) over cosO. The curves are for mass densities of 10 13 , 10 14 , 2 x 10 14 , 3 x 10 14 , 
4 x 10 14 , 5 x 10 14 , 7.5 x 10 14 , and 10 15 g cm 3 . The temperature is 5 MeV, the electron fraction 
is 0.26, the incident electron neutrinos are on their Fermi surfaces, and beta equilibrium has been 
assumed. The 10 13 g cm -3 curve does not include many-body effects, but does incorporate the 
full kinematics. For the other curves, v 2 = 4.5 x 1CT 5 (MeV) -2 , v\ = 1.76 x 10~ 5 (MeV) -2 , and 
an effective mass of 0.75 x m n was assumed, fi is equal to 5.3, 17.4, 26.8, 34.8, 42.1, 48.7, 63.7, 
77.1 MeV for the densities depicted. 



11 



